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1.  Introduction  and  Overview 


The  subject  of  this  report  is  the  characterization  and  interaction 
of  signal  paths  in  very-high-speed  integrated  circuits.  A  long-range 
objective  of  research  on  this  problem  is  optimization  of  signal 
propagation  to  increase  operating  speed  and  minimization  of  both 
internal  and  external  interference.  This  can  be  accomplished  by 
developing  suitable  theoretical  models  for  the  physical  phenomena 
involved.  These  models  are  based  on  circuit  theory  and  linear  system 
theory,  distributed  circuits  and  transmission-line  theory,  and  field- 
theoretic  considerations.  Several  new  and  interesting  results  are 
reported  here  and  many  new  useful  formulas  given. 

First,  the  signal  path  is  modelled  as  a  distributed  circuit  and 
the  propagating  waves  of  voltage,  current,  and  power  flow  (as  a  function 
of  time  and  position  along  the  path)  are  evaluated  analytically  for 
increasingly  complicated,  but  more  accurate,  path  models.  Then  the 
problem  of  matching  the  input  voltage  pulse  shape  to  the  path  or, 
conversely,  the  path  impedance  to  the  pulse  shape  is  considered.  This 
leads  to  the  new  concept  of  nominal  resistance  and  a  simple  way  to 
predict  peak  possible  path  currents.  Next,  the  tangential  components 
of  the  electromagnetic  field  at  the  surface  of  an  integrated  circuit 
(that  couple  to  and  excite  the  microcircuits)  are  derived  for  an 
incident  electromagnetic  wave  or  pulse.  These  fields  as  well  as 
reflected  and  transmitted  wave  power  depend  in  detail  on  the  complex 
reflection  coefficient  which  also  is  given  in  both  exact  and  approximate 
forms.  Then,  the  effects  of  signal  path  electrical  parameters  in 
crosstalk  and  coupling  is  investigated  and  the  coupled  responses  (or 
cross  impulse  and  step  responses)  are  determined  in  analytical  form 
based  on  a  circuit-theoretic  model.  Finally,  the  three-layer  geometry 
of  an  integrated  circuit  excited  by  an  infinitesimal  oscillating 
current  moment  at  the  substrate  surface  is  analyzed  to  determine 
approximate  analytical  expressions  for  the  complete  electromagnetic 
field  at  other  points  on  the  surface.  Aside  from  the  quasistatic 
near  fields,  further  coupling  and  interference  can  occur  due  to  the 
radiation  of  both  a  surface  wave  in  the  substrate  and  a  lateral  wave 
above  the  substrate. 


2.  Summary  of  Research  Results 

Each  of  the  following  six  major  research  results  is  discussed  in 

more  detail  in  a  separate  section  of  this  report. 

•  Complete  analytical  expressions  for  propagating  waves  of 
voltage,  current,  and  power  in  the  time  domain  for  three 
different  transmission-line  models  of  integrated  circuit  signal 
paths. 

•  Determination  of  finite-energy  and  finite-power  voltage  pulse 
shapes  that  are  matched  to  the  signal  path  and/or  load  to  produce 
maximum  current  and  instantaneous  power. 

•  Development  of  the  concept  of  the  nominal  resistance  of  a  complex 
impedance  which  can  be  used  to  predict  maximum  possible  path/load 
currents. 

•  Derivation  of  the  complex  reflection  coefficient  and  tangential 
electromagnetic  fields  at  the  surface  of  an  integrated  circuit 
due  to  an  incident  electromagnetic  wave  or  pulse. 

•  Analysis  of  inductive  and  capacitive  coupling  mechanisms  in 
integrated  circuits  to  determine  coupled  or  cross-impulse  and 
cross  step  responses  and  their  dependence  on  lumped  parameters 
of  the  signal  paths. 

•  Theoretical  study  of  the  complete  electromagnetic  field  produced 
by  an  infinitesimal  oscillating  current  moment  located  on  the 
dielectric  substrate  of  an  integrated  circuit  including  discovery 
of  a  propagating  surface  wave  (in  the  substrate)  and  a  lateral 
wave  (above  the  substrate). 


Technical  Papers: 

"Nominal  Resistance  of  a  Complex  Impedance,"  Electronics  Letters, 
vol.  22,  no.  13,  pp.  703-704,  June  1986. 

"Peak  Transient  Current  and  Power  into  a  Complex  Impedances 
Theory  and  Experiment,"  (in  review),  31  pages,  12  figures,  1  table. 
May  1986. 

"Cross  Impulse  and  Step  Response  Models  for  Capacitive  and 
Inductive  Coupling  in  Very-High-Speed  Integrated  Circuits," 

(in  preparation),  June  1986. 

"The  Propagation  of  Signals  Along  a  Three-Layered  Region,"  (in 
preparation),  August  1986. 
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3.  Distributed-Circuit  Signal  Path  Characteristics 

In  an  integrated  circuit  signal  path  the  circuit  parameters 
(resistance,  capacitance,  inductance,  and  conductance)  are  distributed 
on  a  per-unit-length  basis.  Therefore,  conventional  circuit  analysis, 
which  concentrates  or  ''lumps”  these  parameters  into  circuit  elements, 
is  not  an  accurate  approximation  unless  the  frequency  of  operation  is 
very  low.  If  the  actual  waveshape  of  the  voltage  or  current  that 
propagates  along  the  path  is  required  then  a  distributed  circuit  (or 
transmission-line)  model  of  the  signal  path  is  a  better  approximation. 

The  purpose  of  this  section  is  to  present  analytical  formulas 
derived  for  voltage,  current,  and  power  (energy/time)  as  a  function  of 
time  and  position  on  the  signal  path  for  three  different  distributed 
circuit  models.  The  response  of  these  different  models  to  a  unit  step 
of  input  voltage  is  given.  Path  response  to  an  arbitrary  logic 
sequence  of  voltage  step  transitions  can  be  computed  by  superposing 
step  responses. 

Figure  3.1  is  a  schematic  representation  of  an  integrated  circuit 
signal  path  excited  at  time  t  =  0  by  a  unit  step  function  of  voltage 
(a  transition  at  t  *  0  from  zero  to  one  volt).  The  problem  is  to  find 
the  path  voltage  v(x,  t)  and  current  i(x,  t)  at  a  distance  x  from  the 
source  when  the  distributed  parameters  R,  C,  G,  and  L,  the  series 
resistance,  shunt  capacitance,  shunt  conductance,  and  series  inductance 
(all  per  unit  length)  are  specified. 


Case  1.  G  =  L  =  0.  Here  it  is  assumed  that  series  inductance  and 


shunt  conductance  are  negligible.  (In  an  integrated-circuit  signal 
path  model  these  parameters  are  very  small.)  This  is  called  the  RC 
line  or  Thomson  cable.  The  voltage  v(x,  t)  in  response  to  a  unit  step 
voltage  input  is. 


v(x,  t)  =  erfc 


2 

x  RC 


volts , 


(3.1) 


4t 


where  erfc(y)  is  the  complementary  error  function  of  y.  Rational 
function  approximations  of  the  error  function  and  its  complement  are 
available  in  Handbook  of  Mathematical  Functions  (Abramowitz  and  Stegun, 
eds.)  #7.1.27  and  #7.1.28,  page  299,  The  current  that  flows  in  response 
to  a  unit  step  voltage  input  is, 


i (x,  t)  = 


-x2RC/(4t) 
e  amp. 


(3.2) 


2 

In  (3.1)  and  (3.2)  the  quantity  x  RC  is  the  product  of  total  path 

resistance  and  capacitance  of  the  path  of  length  x.  The  voltage  rises 

2 

to  50%  of  its  final  value  in  approximately  x  RC  seconds;  the  current 

2 

maximum  occurs  at  t  =  0.5  x  RC  and,  the  instantaneous  power  (product 

2 

of  voltage  and  current)  is  a  maximum  when  t  =  1.6  x  RC.  Thus  the 
.  2 

quantity  x  RC  is  like  a  time  constant  that  gives  an  order-of-magnitude 
estimate  of  the  "speed"  of  the  RC  line.  Stated  in  different  terms, 
"slowness"  of  the  signal  path  is  proportional  to  path  length  squared 
and  the  RC  product. 

Case  2.  L  =  0.  If  the  shunt  conductance  G  of  the  signal  path  is 
significant  (due  to  losses  in  the  dielectric  substrate)  then  it  is 
necessary  to  include  a  non-zero  G  in  the  analysis.  This  yields  the 
so-called  "leaky  cable."  There  are  two  loss  mechanisms,  the  resistance 
of  the  conducting  signal  path  itself  and  the  conductance  of 
insulating  substrate  medium.  The  results  of  an  analysis  for  the 
voltage  and  current  in  response  to  a  unit  step  voltage  input 
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I 


and, 


i(x,  t) 


-x2RC/(4t)  , 
e  + 


"\ 

e-x2EC/(4-Oe-Gr/cdT 


(3.4) 

Comparison  of  (3.3)  and  (3.4)  with  (3.1)  and  (3.2)  indicates  that 

there  are  two  effects  of  the  "leakage"  caused  by  shunt  conductance  G. 

First,  the  current  and  voltage  waves  are  attenuated  exponentially  in 

~Gt  /C 

time  by  the  factor  e  .  Secondly,  a  component  is  added  to  each 

response  consisting  of  a  progressive  integral  of  the  attenuated  wave. 

The  major  effect  of  G  is  attenuation  of  signal  strength  due  to  loss. 

A  secondary  effect  is  waveshape  change  of  the  step  responses.  Both  of 

these  effects  can  be  minimized  by  reducing  the  G/C  ratio  or  increasing 

path  capacitance  per  unit  length. 

In  both  Cases  1  and  2  the  absence  of  inductance  L  eliminates  the 

usual  propagation  delay  of  x  <J  LC  seconds  normally  associated  with 

transmission  lines.  Each  of  the  previous  cases  represents  overdamped 

propagation  because  the  ratio  of  attenuation  constant  to  phase  constant 

is  equal  to  or  greater  than  one  so  that  waves  damp  out  in  less  than 

one  wavelength.  Although  they  are  more  complicated,  (3.3)  and  (3.4) 

are  probably  a  better  approximation  to  the  physical  situation  and  could 

be  used,  for  example,  to  "fit"  experimentally  measured  data. 

There  are  two  factors  affecting  overall  "speed"  of  the  signal  path. 

The  first  is  the  slowness  and  attenuation  of  a  step  transition  due  to 

characteristic  times  associated  with  resistance/capacitance  combinations 
2 

(e.g.,  x  RC  and  C/G  seconds).  The  second  factor  is  the  introduction  of 
an  overall  propagation  delay  equal  to  x  f LC  seconds  when  series 
inductance  L  ^  0.  While  the  series  inductance  is  normally  very  small, 
it  does  affect  voltage  and  current  waveshapes,  too.  These  considerations 
lead  to  the  third,  and  most  general,  case  in  which  all  parameters 


R,  C,  G,  and  L  are  non-zero. 

Case  3.  All  parameters  non-zero.  The  voltage  and  current  in 
response  to  a  unit  voltage  step  excitation  of  the  general  line  is 
more  complicated  than  in  the  previous  two  cases  but  still  computationally 
tractable.  Owing  to  a  propagation  delay  of  x  /LC  seconds,  both 
responses  are  zero  until  t  £  x  /LC  when  they  have  the  general  form. 


1/R  G\ 

- —  +  -  x  / LC" 

.  2  L  C/ 

v(x,  t)  =  e 


(3.5) 


1/R  G 


/  R  G\ 


J ;  21 +  m 


1/R  G 


C/1  1 


2  2 
X  -  x  LC 


- --  ♦  -It 


i(x,  t)  =  f  L  e  2'L  Iol2'L  CHt2  -  xLC 


1/R  G 


(3.6) 


1/R  G 


+!  5 


1/R  G 


2\L  C /  lx2  -  x2LC  df  i  , 


x  {LC 


for  t  i  x  |/lC.  The  modified  Bessel  functions  of  the  first  kind  I  (y) 

o 

and  Ij(y)  in  (3.5)  and  (3.6)  are  represented  by  rational  approximations 
in  Handbook  of  Mathematical  Functions  (Abramowitz  and  Stegun,  Eds.) 
#9.8.1  through  #9.8.4,  page  378. 

Two  combinations  of  the  four  line  parameters  appear  in  (3.5)  and 
(3.6).  The  combination  (R/L  +  G/C)  is  an  attenuation  or  damping 
term  while  the  combination  (R/L  -  G/C)  is  proportional  to  the  amount 
of  signal  linear  distortion  the  line  causes  (in  the  distortion-free 
line  R/L  «  G/C,  so  that  the  distortion  terms  are  zero). 
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In  summary,  (3.1)  through  (3.6)  are  complete  analytical 
expressions  for  the  voltage  wave  v(x,  t)  and  current  wave  i(x,  t) 
that  propagate  along  a  transmission-line  signal  path  model  when 
excited  by  a  unit  step  function  of  voltage.  The  corresponding  power 
flow  p(x,  t)  is  the  product  of  these  two  waves  and  equals  the  energy/time 
at  the  position  x  and  at  time  t.  Equations  (3.1)  and  (3.2),  for  the 
RC  line,  are  valid  when  both  the  series  inductance  L  and  shunt  con¬ 
ductance  G  of  the  path  can  be  neglected.  When  substrate  losses  are 
significant,  the  shunt  conductance  G  must  be  included.  This  leads  to 
(3.3)  and  (3.4)  in  which  the  current  and  voltage  waves  of  the  RC  line 
are  attenuated  and  distorted  due  to  the  G/C  ratio.  If  finite  series 
inductance  is  incorporated  then  an  overall  propagation  delay  xvLC~ 
seconds  is  introduced  and  both  the  voltage  and  current  waves  consist 
of  distorted  and  non-distorted  parts.  These  last  results  are  more 
complicated  because  they  involve  integrals  of  Bessel  functions. 

The  three  different  transmission-line  or  distributed  circuit  models 
of  integrated  circuit  signal  paths  presented  here  are  useful  because 
they  give  the  voltage  and  current  waves  in  the  time-domain  (which  can 
be  compared  directly  to  experimental  measurements).  Using  these  models 
the  engineering  tradeoffs  between  pulse  strength,  shape,  propagation 
characteristics  and  ohmic  losses  can  be  computed  for  arbitrary  choices 
of  path  parameters. 

In  practice,  the  actual  signal  path  is  finite  and  terminated 
with  some  load.  This  gives  rise  to  reflected  waves  of  voltage  and 
current.  These  can  be  computed  knowing  the  incident  wave  and  load. 

A  more  rigorous  and  comprehensive  analytical  approach  to  this 
general  problem  would  use  a  field-theoretic  analysis  to  determine  all 
possible  modes  of  wave  propagation  including  a  determination  of  the 
complex  wavenumber  kc  for  the  current  that  propagates  along  the 
conducting  signal  path.  An  analysis  of  this  sort  would  reveal  the 
extent  to  which  distributed-circuit  transmission-line  models  accurately 
represent  dominant  wave  phenomena  on  integrated  circuit  signal  paths. 

An  alternative  approach  is  to  compare  predictions  from  transmission-line 
models  with  experimental  measurements  on  actual  signal  paths. 
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4.  Nominal  Resistance  and  Matching 

Two  important  issues  to  consider  when  devising  an  optimum  signal 
pulse  launch-propagate-receive  system  are:  (1)  how  to  match  a  voltage 
pulse  shape  to  the  input  impedance  of  the  loaded  signal  path  or  to 
the  load  itself  in  order  to  maximize  the  input  current  and  power,  and 
(2)  given  the  voltage  pulse  shape,  how  is  the  optimum  signal  path  or 
load  impedance  determined.  Two  classes  of  voltage  signals  were  con¬ 
sidered.  First,  where  the  total  energy  of  the  input  voltage  is  fixed 
and  its  shape  is  to  be  determined.  Second,  where  the  maximum  amplitude 
of  the  voltage  signal  is  limited  and  its  transitions  from  maximum  to 
minimum  value  are  to  be  determined.  A  secondary  consideration  is 
whether  operation  is  bandlimited.  Both  of  these  problems  have  been 
solved  and  simple  formulas  for  the  voltage  signals  and  their  relation 
to  load  impedance  are  available. 

For  each  class  of  input  signal  (i.e.,  finite  energy  or  finite  power) 
a  new  quantity  called  the  nominal  resistance  of  a  complex  impedance 
was  defined  as  the  ratio  of  peak  applied  voltage  to  peak  current  under 
matched  conditions.  Nominal  resistance  is  a  function  of  the  electrical 
properties  of  the  signal  path  and/or  load,  namely,  the  inverse  Fourier 
transform  of  the  reciprocal  of  the  impedance.  Thus,  knowing  the 
nominal  resistance  and  the  class  of  input  voltage  signal  the  maximum 
(or  peak)  path  or  load  current  can  be  predicted. 

A  summary  of  the  mathematical  details  and  some  supporting  experimental 
data  are  included  in  the  reprint  (see  next  page)  "Nominal  Resistance 
of  a  Complex  Impedance,"  Electronics  Letters,  vol.  22,  no.  13, 
pp.  703-704,  June  1986.  A  more  comprehensive  paper  on  this  general 
subject  entitled,  "Peak  Transient  Current  and  Power  into  a  Complex 
Impedance:  Theory  and  Experiment"  (31  pages,  12  figures,  1  table)  is 
currently  in  review. 
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NOMINAL  RESISTANCE  OF  A  COMPLEX 
IMPEDANCE 

Indexing  terms:  Circuit  theory  and  design.  Time-domain 
analysis.  Impedance 

The  theoretical  relationships  between  broadband  electrical 
input  impedance  and  maximum  instantaneous  input  current 
are  given  for  classes  of  input  voltage  signals  with  or  without 
bandwidth  limitations  A  nominal  resistance  associated  with 
a  complex  impedance  is  defined,  and  this  quantity  can  be 
used  in  circuit  analysis  or  design  to  predict  maximum  pos¬ 
sible  transient  current.  Results  from  several  experimental 
measurements  are  presented. 


When  a  complex  load  impedance  Z(«)  is  driven  by  an  arbi¬ 
trary  voltage  waveform  t> (/),  the  maximum  value  i.,,  of  the 
transient  current  i(t)  that  flows  into  the  load  cannot  be  pre¬ 
dicted  without  carrying  out  a  detailed  mathematical  analysis 
to  obtain  the  transient  response  in  analytical  form.1  This  is 
impractical  w  hen  tit)  is  a  member  of  a  class  of  unpredictable 
waveforms  such  as.  for  example,  a  logic  sequence  of  step  tran¬ 
sitions.  a  noise  process  or  finite-energy  analogue  signals  corre¬ 
sponding  to  segments  of  speech  or  music.  Prior  knowledge  of 
is  essential,  however,  when  specifying  the  output  current 
and  power  capabilities  of  the  circuit  that  produces  the  voltage 
tHf )  to  drive  the  load  impedance  Z(a>).  Three  specific,  practical 
examples  are  (il  estimating  interstage  peak  currents  in  ampli¬ 
fiers.  (ii)  predicting  logic-line-driver  peak  currents  in  inte¬ 
grated  circuits,  and  (iii)  determining  audio-power-amplifier 
peak  output  current  required  to  drive  a  given  loudspeaker 
system. 

In  this  letter  a  new  quantity  is  introduced,  called  the 
nominal  resistance  associated  with  Z(ai),  defined  as  the  ratio 
of  peak  applied  voltage  to  peak  current.  Nominal  resistance 
depends  not  only  on  Z(iu|  but  also  on  the  class  of  input  signal 
(finite  energy  or  finite  power!  and  the  operating  bandwidth. 
Knowing  the  nominal  resistance,  the  maximum  possible  peak 
current  can  be  predicted  for  design  or  analysis  purposes. 

The  current  i(/>  that  flows  into  a  complex  impedance  Z(ui( 
in  response  to  a  voltage  il/)  is  given  by  the  temporal  convolu¬ 
tion 


x 

=  j"  iftlili  -  x) 


.  i  J 


do 


which  is  the  inverse  Fourier  transform  of  the  admittance 
y(cu)  =  I  Z(u>).  The  quantity  ylr)  can  also  be  interpreted  as  the 
current  in  response  to  a  unit  impulse  of  voltage. 

Two  separate  problems  are  considered  here.  First,  maximise 
the  current  i(r|  when  the  total  energy  £  of  ilf)  is  finite: 


x 

-J' 


2(/|  di  <  y. 


that  is.  find  the  shape  of  ill).  Secondly,  maximise  i(fl  when  the 
amplitude  of  tl/)  is  restricted  by 

0  <  MO  |  £  i0  <  x  (4) 

where  r0  is  a  constant.  It  is  assumed  that  eqn  2  exists  and 
that  if  ylr)  contains  singularity  functions  they  are  of  order  no 
higher  than  a  delta  function  A/I.  Thus,  for  example,  Z(ui)  may 
not  contain  a  capacitor  across  its  input  terminals. 

The  solution  to  the  first  problem,  from  the  matched-filter 
principle,2  is 

t.OI  =  M-z)  (5) 


where  the  constant  k  is  determined  from  eqn.  3.  The  peak 
current  squared,  in  response  to  eqn.  5,  is2 

«  a 

i2(0)=£  |  y2(t)dt  =  ~  |  |  V(w)  |2  dai  (6) 


with  £  specified  by  eqn.  3.  The  ratio  of  peak  voltage  to  peak 
current  or  nominal  resistance  for  the  matched-filter  input 
voltage  eqn.  5,  found  from  combining  eqns.  3,  5  and  6,  is 

max  [></)] 

Rm  =  Wl  y2(/)  di  ,7' 

For  the  second  problem,  the  convolution  eqn.  1  is  maximised 
when  the  input  voltage  is 

rjf)  =  i'0  sign  [><-/)]  (8) 

The  peak  current  in  response  to  eqn  8  is 


x 

=  t0  J  Id'll  di 


The  nominal  resistance,  corresponding  to  eqns.  8  and  9,  is 


•  f- ,  Id'll  dl  llu' 

If  operation  is  limited  to  the  finite  bandwidth  a  rad  s.  then 
substitution  of  r„(r)  for  ylr)  in  eqns.  7  and  10  is  necessary, 
where  »•„(/)  is  computed  from  eqn.  2  using  finite  integration 
limits  ±o. 

Fig.  1  illustrates  impedance  magnitude  against  frequency 
for  a  14  35  fi  carbon  resistor  R  and  five  different  complex 
impedances  numbered  1-5.  These  data  were  measured  digi¬ 
tally  and  represent  discrete  values  every  48  Hz  throughout  a 
20  kHz  bandwidth  Table  1  lists  both  the  DC  value  (u  =  0) 
and  minimum  value  (for  to  >  0)  of  the  impedance  magnitudes 
in  Fig.  I.  The  quantity  yin  in  eqn.  2  was  estimated  for  each 


Fig.  I  Input  impedance  magnitude  of  Id  33  O  resistor  R  and  hie  differ¬ 
ent  complex  impedances 

impedance  as  the  current  in  response  to  a  20  /is-wide  rec¬ 
tangular  voltage  pulse.  Using  1024  samples  (50  kHz  sampling 
frequency!  of  each  yin.  based  on  an  average  of  2000  trials  to 
improve  signal  noise  ratio,  nominal  resistances  were  evaluated 
numerically  by  approximating  the  integrals  in  eqns  7  and  10 
with  finite  sums.  The  resulting  numerical  values  of  R ,  and  Rm 
are  also  given  in  Table  1. 

Examination  of  the  data  in  Table  1  will  reveal  thai  Rm.  the 
nominal  resistance  corresponding  to  a  matched-filtered.  fimie- 
energy  input  voltage,  was  not  found  to  be  less  than  the 
minimum  impedance  magnitude,  but  R,.  the  nominal  resist¬ 
ance  for  a  voltage-limited  finite-power  input,  can  be  approx¬ 
imately  half  of  the  minimum  impedance  magnitude. 

The  accuracy  of  a  nominal  resistance  computation  depends 
on  how  well  the  quantity  yir)  in  eqn.  2  is  estimated  experimen- 
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Table  1 


Impedance 

2(0)  —  Rpc 

IZMLi, 
to  >  0 

R. 

n 

o 

a 

a 

R 

14  35 

14  35 

1435 

14  35 

1 

8  84 

938 

10-86 

6  12 

2 

6  18 

4-67 

689 

3  12 

3 

705 

295 

3-84 

1  51 

4 

3-57 

3  28 

6  18 

1  86 

5 

20-66 

8  38 

11-97 

467 

tally.  Comparison  of  the  inverse  FFT  of  the  sampled  values  of 
><r)  derived  from  the  20  /js  pulse  measurements  indicated  that 
the  admittance  spectral  magnitude  |  V'(u>)|  was  within  ±  I  dB 
of  the  true  admittance  magnitude  measured  independently.  In 
this  connection,  it  also  is  possible  to  estimate  y(r)  using  differ¬ 


ent  voltage  pulse  shapes  or,  indirectly,  as  the  inverse  FFT  of 
the  measured  complex  admittance. 
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5.  Peak  Path  Currents  and  Interference 

The  concept  of  nominal  resistance  introduced  in  the  previous  section 
has  further  application  in  the  study  of  crosstalk,  inductive  coupling, 
and  interference  between  signal  paths.  For  a  given  signal  path 
and/or  load  there  exists  a  specific  logic  sequence  that  will  excite 
extremely  large  peak  instantaneous  currents  (as  well  as  large  rates  of 
change  of  current).  Such  currents  can  potentially  interfere  with 
other  circuits  through  the  mechanism  of  inductive  coupling.  If  the 
input  impedance  of  the  path  and/or  load  is  Z(w)  then  the  logic  sequence 
producing  the  maximum  possible  value  of  instantaneous  current  is  produced 
by  the  voltage  v(t)  =  vosign[y(-t) ]  where  vq  is  the  maximum  voltage 
of  the  logic  sequence  and  y(t)  is  the  inverse  Fourier  transform  of 
Y(co)  =  l/Z(«o).  The  peak  current  and  nominal  resistance  formula  are 
given  in  equations  (9)  and  (10),  respectively,  of  the  reprint  in  the 
previous  section. 

When,  in  analysis  or  design  of  integrated  circuit  signal  paths 
and  loads,  the  question  of  maximum  possible  instantaneous  current  arises, 
the  answer  is  simple  and  direct.  The  maximum  possible  transient  current 
Snax  *n  resP°nse  to  anY  logic  sequence  is  given  by  a  formula  similar  in 
form  to  Ohm's  Law,  namely,  i  *  v  /R  .  In  this  equation,  R  is  the 
nominal  resistance  of  the  path  and  v^  the  maximum  voltage  of  the  logic 
sequence.  Thus  nominal  resistance  is  an  important  new  parameter  that 
characterizes  a  signal  path  and/or  load.  The  smaller  Rr  is  the 
larger  the  maximum  transient  current  can  be.  Nominal  resistance  can  be 
c; lculated  from  y(t),  the  current  in  response  to  a  unit  impulse 
excitation  of  the  signal  path,  using  equation  (10)  in  the  previous 
section.  The  currents  i(x,  t)  in  (3.2),  (3.4),  and  (3.6)  of  section  3 
of  this  report  can  be  used  to  compute  Rn  for  transmission  line  models 
of  signal  paths  if  they  are  first  differentiated  with  respect  to  time 
(because  they  are  the  currents  in  response  to  a  unit  step  rather  than 
impulse  of  voltage). 
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6.  Incident  Electromagnetic  Waves  or  Pulses 

In  terms  of  high-frequency  electromagnetic  fields  an  integrated 
circuit  is  a  three-layer  problem  consisting  of  air  above,  a  lossy 
dielectric  substrate,  and  an  imperfectly  conducting  metal  ground  plane 
underneath.  When  this  circuit  is  illuminated  from  above  by  an  incident 
electromagnetic  plane  wave  (or  an  electromagnetic  pulse  consisting  of  a 
superposition  of  plane  waves  of  different  frequencies)  part  of  each  wave 
is  reflected  away  and  the  remainder  transmitted  (refracted)  into  the 
dielectric  substrate  and  conducting  layers. 

An  important  quantity,  which  completely  characterizes  the  reflection 
of  incident  electromagnetic  waves  at  the  air-dielectric  interface,  is 
the  complex  reflection  coefficient  f^.  The  squared  magnitude  of  the 
reflection  coefficient  0  ^  [f^^S  1  is  a  direct  measure  of  the  ratio 
of  reflected  power  to  incident  power.  The  reflection  coefficient  f 
depends  on  the  electromagnetic  properties  of  each  of  the  three  layers 
(i.e.,  permittivity,  permeability,  and  conductivity),  the  incident  wave’s 
frequency,  and  the  orientation  of  the  incident  wave's  electric  and 

magnetic  vectors  with  respect  to  the  plane  of  incidence. 

2 

For  a  given  three-layer  geometry,  ( f r  |  indicates  what  proportion  of 
the  incident  wave's  power  is  reflected  (harmlessly)  away  while 
(potentially  damaging)  wave  power  that  is  absorbed  by  the  integrated 
circuit  is  1  -  [f^  J  .  Therefore,  to  minimize  possible  damage  to  an 
integrated  circuit  caused  by  an  incident  electromagnetic  wave  of  frequency 
w,  the  reflection  coefficient  at  that  frequency  should  be  |fr  |~1 
for  all  likely  angles  of  incidence.  When  dealing  with  an  incident, 
high-energy  electromagnetic  pulse,  the  potentially  damaging  effects 
of  electromagnetic  radiation  are  minimized  when  |fr  |  ~  1  for  those 
frequency  ranges  throughout  which  the  spectrum  of  the  pulse  contains 
significant  energy. 

The  reflection  coefficient  also  provides  insight  into  the  problem 
of  shielding  or  radiation  hardening  of  integrated  circuits.  Because  the 
three-layer  geometry  and  its  electromagnetic  properties  as  well  as  the 
frequency  range  of  incoming  electromagnetic  energy  (usually)  are  all 
given  quantities,  the  frequency  ranges  where  the  reflection  coefficient 


|  f r  I 0,  or  small  compared  with  unity  are  vulnerable  and,  consequently, 
do  need  some  external  shielding. 

The  fact  that  f  is  a  complex  function  (i.e.,  having  both 
magnitude  and  phase)  has  further  technical  implications.  Certain 
incident  pulse  shapes  may  produce  extremely  high  momentary  field 
strengths  in  the  integrated  circuit  if  they  are,  for  example,  conjugate 
matched  to  the  complex  transmission  coefficient  (which  is  related  to 
the  reflection  coefficient).  This  situation  is  similar  to  matched 
filtering.  Due  to  phase  linearization,  waves  superpose  at  one  instant  in 
time  rather  than  being  dispersed  in  time.  Thus,  information  provided  by 
the  complex  transmission  coefficient  is  useful  in  two  ways.  First,  to 
shield  against  such  pulses  and,  secondly,  to  design  electromagnetic 
pulse  shapes  that  are,  for  a  given  energy,  potentially  most  damaging  to 
integrated  circuits. 

A  further  use  of  the  complex  reflection  coefficient  is  in  the 
determination  of  the  tangential  components  of  the  electric  and  magnetic 
field  at  the  interface  between  air  and  dielectric  surface.  These  fields 
can  excite  (couple  to)  the  circuits  on  the  dielectric  surface. 

Figure  6.1  shows  the  three-layer  geometry  consisting  of  air 
(Region  0),  dielectric  (Region  1),  and  conductor  (Region  2).  Consider  a 
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Fig.  6.1  Plane  wave  incident  on  boundary  z  =  0  of  dielectric 
Region  1. 
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plane  wave  incident  from  the  air  (Region  0,  z  ^  0)  on  the  plane  boundary 
z  =  0  between  regions  0  and  1  as  shown  in  Fig.  6.1.  Region  1  has 
uniform  thickness  1  and  is  bounded  at  z  =  i  by  Region  2  (z2J?  ). 

The  complex  wave  number  in  each  region  is. 


1/2  , 


kj  =  p.  +  '*j  =a,[/0(€j  +  ,  j  =  0,  1,  2, 


(6.1) 


where  €  and  tr  are  the  real,  effective  permittivity  and  conductivity, 


respectively,  and  the  permeability/^  is,  by  assumption,  that  of  air. 


The  angles  9 Q,  0^,  and  6^  in  Fig.  6.1  are  related  by  Snell's  law. 

For  an  electromagnetic  field  of  the  electric  type  (i.e.,  with  the 
electric  field  in  the  plane  of  incidence),  the  result  of  a  lengthy 


derivation  for  the  complex  reflection  coefficient  f  is, 
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Kith  the  notation. 
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k  =  xk  +  zk  ,  k 
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=  ksinG  and  k  =  kcos#. 


(6.3) 


the  reflection  coefficient  can  be  expressed  as. 
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This  is  the  general  formula  for  the  plane-wave  reflection  coefficient  of 
the  electric  type  for  the  three-layered  region. 

The  general  formula  (6.4)  can  be  simplified  with  the  practically 
valued  inequalities, 
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to  the  form. 
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In  this  approximation  the  first  term  is  the  reflection  coefficient  of 
the  two-layered  region  with  Jl  =  0. 

The  electromagnetic  field  of  electric  type  on  the  air-surface  of 
the  dielectric  is. 
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o 

where  f^  is  given  by  (6.4)  or  (6.6).  The  tangential  components 

B  (x,  0)  and  E  (x,  0)  are  of  primary  interest  in  the  excitation  of 
oy  ox 

circuits  on  the  dielectric  surface. 

The  reflection  coefficient  for  fields  of  the  magnetic  type 
(i.e.,  with  the  magnetic  field  in  the  plane  of  incidence)  for  the 
three-layer  problem  of  Fig.  6.1  is, 
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Kith  the  notation. 


k  =  ksin0,  k  =  kcos©, 
x  z 


(6.11) 


(6.10)  becomes, 
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(6.12) 


This  is  the  general  form  for  the  plane-wave  reflection  coefficient  of  the 
magnetic  type  for  the  three-layered  region. 

Kith  the  condition. 
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(6.12)  simplifies  to  the  form. 
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In  this  approximation  the  first  term  is  the  reflection  coefficient  of 
the  two-layer  region  obtained  with  £  =  0. 

The  electromagnetic  field  of  the  magnetic  type  on  the  surface  of  the 
dielectric  is, 


ik  x 

E  (x,  0)  =  E  lnC(0,  0)e  °X  [1  +  f  ], 
oy  oy  mr 


(6.15) 
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Boz(x’  0)  =  77  EoyinC(0*  0)elk°XX[1  +  fmr]sin*o’ 


(6.17) 


The  tangential  components  E  (x,  0)  and  Bqx(x,  0)  play  an  important 
role  in  exciting  the  electrical  circuits  on  the  dielectric  surface. 

In  summary,  the  following  new  research  results  are  contained  in 
this  section.  For  a  plane  wave  of  arbitrary  angle  of  incidence  on  the 
dielectric  surface  of  a  three-layer  model  for  an  integrated  circuit, 
the  complex  reflection  coefficient  is  given  in  both  exact  and  approxi¬ 
mate  analytical  form  for  the  electric  field  of  the  wave  in  the  plane 
of  incidence  [(6.4)  and  (6.6)]  and  for  the  magnetic  field  of  the  wave 
in  the  plane  of  incidence  [(6.12)  and  (6.14)].  The  tangential  components 
of  the  electromagnetic  field  capable  of  exciting  circuits  at  the  air- 
dielectric  interface  are  given  in  exact  analytical  form  [(6.7),  (6.8), 
(6.9)  and  (6.15),  (6.16),  (6.17)].  These  expressions  are  relatively 
simple  and  suitable  for  numerical  computation.  Because  frequency  go  is 
an  arbitrary  parameter  contained  in  the  complex  wavenumbers  (6.1),  these 
results  also  are  applicable  to  electromagnetic  pulses  composed  of  a 
superposition  of  (or  spectrum  of)  plane  waves.  The  general  formulas 
can  be  used  to  analyze  both  existing  geometries  and  materials  or  proposed 


new  ones 


7.  Circuit-Theoretic  Path  Coupling  Model 


The  objective  of  this  section  is  to  present  analytical  formulas 
and  curves  for  coupling  between  signal  paths  in  integrated  circuits 
based  on  lumped-parameter  circuit  theory  models.  UTiile  this  is  an 
oversimplification  of  the  general  coupling  problem,  the  conclusions 
based  on  the  analysis  are  consistent  with  those  drawn  from  the  more 
complicated  field-theoretic  path  coupling  model  in  section  8  of  this 
report. 

Simple,  capacitively  coupled  RC  circuits  are  analyzed  to  illustrate 
the  time-domain  and  frequency-domain  nature  of  capacitive  coupling. 

Next,  a  more  general  coupled  circuit,  that  includes  both  inductance  and 
resistive  loading  of  the  signal  path,  is  analyzed  to  show  the  effects  of 
capacitive  and  inductive  coupling  separately.  Analytical  expressions 
are  given  for  coupling  gain  and  coupled  (or  cross)  impulse  and  step 
responses.  Several  curves  of  cross  impulse  and  step  responses  are 
presented.  Means  of  reducing  coupling  are  discussed. 

Consider  the  elementary  circuit  in  Fig.  7.1  used  to  model  first-order 

capacitive  coupling  between  signal  paths  on  an  integrated  circuit.  The 

lower  RC  line  is  driven  by  voltage  v^(t)  and  coupled  via  the  mutual 

capacitance  to  the  upper  RC  line  which  is  not  driven  but  has  an 

output  voltage  v^(t)  due  to  mutual  coupling.  The  problem  is  to  determine 

the  voltage  transfer  function  between  v.  and  v  in  the  frequency  domain, 

mm 

then  find  v^  when  v^  is  either  an  impulse  or  step  function  of  voltage. 

Assuming  C  <  <  C,  the  current  through  C  is  small  relative  to 
m  m 

that  through  C  in  the  lower  line  and  this  leads  to  the  voltage  transfer 
function  (in  terms  of  the  complex  frequency  variable  s). 
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C  1 
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V.  (s)  1  +  sRC  /  1  +  sRC  C 
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(7.1) 


The  term  in  front  of  the  brackets  is  the  voltage  across  C  in  the  lower 

line  while  the  term  in  the  brackets  represents  the  voltage  coupled 

through  C  . 

m 

The  magnitude  of  the  spectrum  of  the  transfer  function  (7.1)  is  shown 
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Fig.  7.1.  First-order  capacitive  couplings  circuit,  complex  transfer 

function  and  coupling  spectrum,  coupled  impulse  response  h(t), 
coupled  step  response  a(t). 
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in  Fig.  7.1.  Coupling  is  bandpass  having  a  maximum  at  the  frequency 

Co  =  =  1/RC  where  the  coupling  coefficient  equals  0.25  C^/C.  To  either 

side  of  cu  =  coupling  rolls  off  at  6dB/oct.  The  inverse  transform 

of  (7.1)  yields  the  "cross-impulse  response"  h(t),  that  is, 

h(t)  =  v  (t)  when  v.  (t)  =  5"(t).  When  the  lower  line  is  excited  by  an 
m  in 

impulse  S(t)  the  coupled  voltage  h(t)  has  the  shape  shown  in  Fig.  7.1. 

The  peak  value  of  h(t)  equals  °<C  /C  and  it  occurs  when  t  =  0.  The 

m 

areas  under  the  positive  and  negative  portions  of  the  impulse  response 
are  equal  and  opposite.  This  occurs  because  the  step  response  a(t), 
which  equals  the  integral  of  h(t),  has  a  final  value  of  zero  (corre¬ 
sponding  to  the  zero-frequency  value  of  the  coupling  transfer  function). 

The  maximum  value  of  the  step  response  is  0.368  C  /C  and  occurs  at 

m 

t  =  l/<*  =  RC. 

As  Fig.  7.1  summarizes,  the  maximum  of  the  voltage  coupling 
coefficient,  cross  impulse  and  step  responses  are  all  proportional  to 
C  /C,  the  ratio  of  coupling  capacitance  to  line  capacitance.  This  ratio 
could  be  made  small  by  increasing  the  capacitance  of  the  signal  path, 
or  decreasing  the  mutual  capacitance,  or  both. 

A  more  realistic  model  would  include  some  resistive  loading  at 
the  end  of  each  signal  path  as  well  as  path  inductance  and  inductive 
coupling.  This  increases  both  the  order  of  the  differential  equations 
describing  the  circuits  and  the  number  of  circuit  parameters. 

The  objective  of  the  second  modelling  approach  to  follow  is  to 
determine  the  voltage  coupled  from  one  signal  path  to  another  due  to 
capacitive  and  inductive  coupling.  Each  coupling  mechanism  is  con¬ 
sidered  separately  and  both  frequency-domain  and  time-domain  responses 
are  evaluated  in  analytical  form.  This  lumped-parameter  signal-path 
model  provides  considerably  more  insight  into  the  nature  of  coupling. 

It  quantitatively  describes  the  dependence  of  coupled  voltages  on 

the  parameters  of  the  signal  path,  namely,  path  resistance,  capacitance, 

inductance,  loading  and  mutual  capacitance  and  inductance. 

Figure  7.2  shows  the  second  lumped-parameter  signal-path  model. 

The  lower  and  upper  signal  paths  each  have  path  resistance  R,  inductance 

L,  and  capacitance  C  and  each  path  is  loaded  with  a  resistor  R  . 

o 

Capacitive  coupling  between  the  two  paths  is  taken  into  account  by  the 
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mutual  capacitance  C  and  inductive  coupling  by  mutual  inductance  L  . 

m  m 

Only  five  circuit  parameters  are  needed  to  describe  coupling: 

R,  L,  C,  Rq,  and  for  capacitive  coupling  or  R,  L,  C,  Rq  and  for 
inductive  coupling.  In  Fig.  7.2,  the  lower  signal  path  is  excited  with 


the  arbitrary  voltage  \\n(s),  where  s  =  <T  +  joj  is  the  complex  frequency 
variable,  and  the  voltage  coupled  to  the  (undriven)  upper  signal  path 
is  V  (s). 

m  C 

To  denote  capacitive  coupling  the  notation  V  (s)  =  V  (s)  is  used. 
— t - t m  m 

The  voltage  transfer  function  V  es)/Y.  (s)  can  be  evaluated  simply 

m  in 

when  it  is  assumed  that  the  current  through  C  is  small  compared  to 

m 

that  through  the  parallel  combination  of  C  and  Rq  in  the  lower  path. 
Also,  the  mutual  inductance  is  ignored.  The  result  is, 
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The  voltage  transfer  function  is  bandpass  with  maximum  mid-frequency 
gain  given  by. 
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At  very  high  frequencies  the  gain  rolls  off  at  -12dB/octave  and  has 
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(7.4) 


At  low  frequencies  the  gain  rises  at  +6dB/octave  and  takes  the  form, 
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These  results  show  that,  for  fixed  path  parameters  R,  L,  and  C, 


capacitive  coupling  is  reduced  when  either  the  mutual  capacitance  C 


and/or  path  load  resistance  Rq  _i£  decreased .  When  path  parameters  are 
adjustable,  capacitive  coupling  is  reduced  by  increasing  the  path 
capacitance  C,  keeping  the  mutual  capacitance  fixed,  so  that  the  ratio 


C  /C  decreases, 
m 


Transient  responses  can  be  derived  from  the  general  voltage  transfer 
function  (7.2)  by  inverse  Laplace  transformation.  Two  transient  responses 
are  of  particular  interest:  the  coupled  voltage  when  the  input  is 


an  impulse  h(t)  or  unit  step  u(t).  These  coupled  voltages  (or  cross 


C  C 

impulse  and  step  responses)  are  denoted  by  hm  (t)  and  a^  (t),  respectively 


When  the  cross  impulse  response  is  known  the  coupled  response  to  an 
arbitrary  input  voltage  waveform  can  be  evaluated  by  temporal  convolution. 
The  cross  step  response  yields  the  actual  time-domain  waveform  of  the 


coupled  voltage  due  to  a  typical  step  transition  of  a  logic  sequence. 


The  impulse  response  h^  (t)  is  the  inverse  Laplace  transform  of 


the  voltage  transfer  function  (7.2)  which  can  be  rewritten  as, 


Vm  (S)  C  . 
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where , 


There  are  three  possibilities  for  the  numerical  values  of  the  parameters 
and  !  complex  conjugates,  real  and  equal,  or  real  and  distinct 
(depending  on  whether  the  quantity  under  the  radical  sign  is  negative, 
zero,  or  positive,  respectively).  Here,  the  last  possibility  is  assumed, 
namely, 

(-  +  — -fi  +  — l 

\  L  RC/  LC  V  R 

o  o' 


Physically,  this  assumption  implies  that  the  coupling  mechanism  is 
overdamped  so  that  the  transient  responses  will  consist  of  exponential 
terms  of  the  form  e  or  te  rather  than  exponentially  damped 
sinusoids. 

When  the  voltage  transfer  function  is  expanded  in  a  partial  fraction 
expansion  and  its  inverse  Laplace  transform  taken,  the  cross  impulse 
response  is  found  to  be, 
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where  the  constants  a^  through  d^  are, 
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with«x'1>  and  <^3  given  in  (7.7),  (7.8),  and  (7.9),  respectively. 

Figures  7.3  and  7.4  are  graphs  of  the  cross  impulse  response  (7.11) 
for  several  choices  of  parameters.  The  step  response  equals  the  integral 
with  respect  to  time  of  the  impulse  response  or  the  inverse  Laplace 
transform  of  1/s  times  the  voltage  transfer  function  (7.2).  It  is, 


m  1  -V  "V 

a  (t)  =  —  .  —  |  (a_  +  b_t)e  +  (c0  +  d~t)e 
m  C  LC  2  2  2  2 


(7.16) 


where  the  constants  a^  through  are, 
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(7.18) 
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withex^,  0^2 »  and  ex^  given  in  (7.7),  (7.8),  and  (7.9)  respectively. 

Figures  7.5  and  7,6  are  graphs  of  the  cross  step  response  (7.16) 
for  the  same  choices  of  parameters  as  in  Figs.  7.3  and  7.4.  In  the 
denominator  of  every  term  of  each  transient  response  is  the  constant 
factor. 
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(7.21) 


The  maximum  value  of  the  step  response  can  be  estimated  to  be  of  the 
order. 
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It  is  interesting  to  compare  this  last  result  with  the  previously  calcu¬ 
lated  value  of  mid-band  gain  (7.3)  which  was. 
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To  minimize  capacitive  coupling,  these  last  two  time-domain  and  fre¬ 
quency-domain  expressions  imply  that  the  product. 
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should  be  made  as  small  as  possible.  If  R/L>>  1  / ( R  C)  then  this 

o 

product  simplifies  to, 
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(7.25) 


which  has  the  interesting  interpretation  of  C^/C  times  the  ratio  of  the 
L/R  time  constant  to  the  RC  time  constant.  To  minimize  reduce 
C^/C,  the  ratio  of  mutual  capacitance  to  line  capacitance,  decrease 
the  L/R  time  constant,  and  increase  the  RC  time  constant.  The  product 
P^,  can  also  be  written  as, 


(7.26) 


and  this  suggests  that,  for  both  fixed  path  ohmic  loss  (or  heating) 

and  mutual  capacitance  C^,  path  capacitance  should  be  increased  and  path 

inductance  decreased  to  reduce  capacitive  coupling. 

To  denote  inductive  coupling  the  notation  V  (s)  =  V  L(s)  is  used 
. — ■  *  '  »  ^  m  m 

in  Fig.  7.2.  The  voltage  transfer  function  V  cs)/V.  (s)  can  be 

m  in 

evaluated  by  assuming  that  the  current  I(s)  in  the  lower  path  is 
insignificantly  influenced  by  I'(s),  the  current  induced  in  the  upper 
path  due  to  mutual  inductance  L^.  Also,  the  mutual  capacitance  is 
ignored.  Under  these  assumptions,  the  voltage  transfer  function  is, 
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The  voltage  transfer  function  (7.27)  is  bandpass  with  maximum 
mid-frequency  gain  given  by, 
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At  very  high  frequencies  the  gain  rolls  off  at  -12  dB/octave  and  has 
the  general  form, 
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At  low  frequencies  the  gain  rises  at  +6  dB/octave  and  takes  the  form, 
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These  results  show  that,  for  fixed  path  parameters  R,  L,  and  C, 

inductive  coupling  is  reduced  when  either  the  mutual  inductance  L^ 

and/or  load  resistance  Ro  i^s  decreased.  When  path  parameters  are 

adjustable,  inductive  coupling  is  reduced  by  increasing  the  path 

inductance  L,  keeping  mutual  inductance  fixed,  so  that  the  ratio 

L  /L  decreases, 
m 

Transient  responses  can  be  derived  from  the  general  voltage  transfer 
function  (7.27)  by  inverse  Laplace  transformation.  As  before,  the 
coupled  voltage  when  the  input  is  an  impulse  h(t)  or  unit  step  u(t) 
is  evaluated.  These  coupled  voltages  (or  cross  impulse  and  step 
responses)  are  denoted  by  h^^Xt)  and  a^(t),  respectively.  Knowing  the 
impulse  response,  the  coupled  response  to  an  arbitrary  input  voltage 
waveform  can  be  evaluated  by  temporal  convolution.  The  step  response 
yields  the  actual  time-domain  waveform  of  the  coupled  voltage  due  to  a 
typical  step  transition  of  a  logic  sequence. 

The  impulse  response  h  ^(t)  is  the  inverse  Laplace  transform  of 
the  voltage  transfer  function  (7.17)  which  can  be  rewritten  as, 
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where. 


=  1/(H0C), 


(7.32) 


and and  <=*,  the  same  as  for  capacitive  coupling  and  given  by  (7.7) 
and  (7.8).  Using  a  partial  fraction  expansion  and  performing  termwise 
inverse  Laplace  transformation,  in  the  same  way  as  for  capacitive 
coupling,  the  coupled  impulse  response  due  to  inductive  coupling  has 
the  general  form, 


(7.33) 


where  the  constants  a^  through  d^  are, 
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No  curves  of  the  cross  impulse  response  (7.33)  need  be  plotted  because 
they  are  identical  to  those  in  Figs.  7.3  and  7.4  when  the  substitutions 

"VL-'V0  andV^3  are  made.  To  verify  these  substitutions 
compare  (7.33)  to  (7.11).  Note,  however,  that  the  polarity  of  the 
inductive  cross  impulse  response  is  inverted  relative  to  the  capacitive 
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cross  impulse  response. 

The  coupled  step  response  due  to  inductive  coupling  has  the  form 
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where  the  constants  through  d^  are, 
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As  with  the  cross  impulse  responses,  the  cross  step  responses  for 

inductive  coupling  can  be  seen  in  Figs.  7.5  and  7.6  with  the  substitutions 

-L  / L— *C  /C  and  ex’.  — ►  oc .  Again,  polarity  is  reversed, 
mm  **  J 

Similar  to  capacitive  coupling,  the  denominator  of  each  term  in 
(7.23)  has  the  constant  factor, 
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This  implies  that  the  maximum  value  of  the  step  response  can  be  estimated 
to  be  of  the  order, 
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This  last  result  can  be  compared  with  the  previously  calculated  midband 
gain  (7.28)  which  was. 
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To  minimize  inductive  coupling,  these  last  two  time-domain  and 
frequency-domain  expressions  imply  that  the  product. 
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generally  should  be  made  as  small  as  possible.  If  R/L  >>  1/(RqC)  then 
this  product  simplifies  to, 
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which  is  the  product  of  L  /L  and  the  ratio  of  the  L/R  time  constant 

m 

to  the  RC  time  constant.  The  product  PT  can  also  be  written  as, 
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P.  =  - .  (7.48) 
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and  this  suggests  that,  for  both  fixed  path  ohmic  loss  (or  heating)  and 
mutual  inductance  Lm>  path  capacitance  should  be  increased  to  reduce 
inductive  coupling. 

The  conclusion  based  on  the  foregoing  analysis  of  the  coupling 
model  is  that  there  are  several  options  available  to  reduce  coupling 
and  interference  in  integrated  circuit  signal  paths:  (1)  increase 
path  capacitance;  (2)  decrease  mutual  capacitance  and  mutual  inductance; 

(3)  increase  path  loading;  (4)  if  ohmic  losses  and  heating  are  not 
critical,  increase  path  resistance.  Option  (1),  increasing  path 
capacitance, appears  most  promising.  This  can  be  accomplished  by  increasing 
the  relative  dielectric  of  the  substrate  and/or  reducing  the  thickness 
of  the  substrate  layer. 
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8.  Field-Theoretic  Path  Coupling  Model 

An  important  type  of  electric  circuit  consists  of  conducting  paths 
on  a  dielectric  sheet  that  is  bounded  on  one  side  by  air,  and  on  the 
other  side  by  a  conducting  plane.  This  is  the  case  for  the  majority  of 
chip  area  in  an  integrated  circuit  because  of  the  many  interconnections. 
The  interaction  of  circuit  elements  on  such  a  structure  includes  not 
only  the  conducting  paths  but  also  the  coupling  due  to  electromagnetic 
waves  that  are  generated  and  travel  in  and  along  the  surface  of  the 
dielectric.  This  is  particularly  important  where  extremely  high 
frequencies  are  present  due  to  very  high  clocking  rates.  It  is  well 
known  that  a  dielectric-coated  metal  plane  supports  a  surface  wave 
that  propagates  in  the  dielectric.  However,  such  surface  waves  are 
usually  derived  as  possible  solutions  of  Maxwell's  equations  when  the 
conditions  of  a  surface  wave  are  imposed  (viz.,  that  the  surface  wave 
be  exponentially  attenuated  in  the  air  in  the  direction  perpendicular 
to  the  dielectric).  It  is  not  obvious  that  such  a  wave  is,  in  fact, 
generated  by  alternating  currents  located  on  the  air-dielectric  boundary. 
A  first  step  in  the  determination  of  the  fields  generated  by  such 
currents  is  a  complete  analysis  of  the  electromagnetic  field  of  a 
unit  infinitesimal  horizontal  dipole  (i.e.,  an  infinitely  small  con¬ 
ductor  with  an  electric  current  moment  Ihe  =  1  ampere-meter)  located 
on  the  air-dielectric  surface  as  shown  in  Fig.  8.1. 


(a)  (b) 


Fig.  8.1(a)  Side  view  of  unit  electric  dipole  on  plane  boundary 
z  =  0  between  air  and  a  dielectric;  (b)  Top  view  of  unit  electric 
dipole  on  dielectric  surface  and  observation  point  P(^~,  0)  where 

the  o  and  <$>  components  of  the  electric  field  vector  are  to  be 
evaluated. 


-41- 


Once  this  is  known,  the  fields  generated  by  conductors  of  finite 
length  can  be  derived.  These  include  single  conductors,  parallel 
conductors  forming  a  transmission  line,  and  various  circuits  with 
small  lumped  elements. 

The  electromagnetic  field  generated  by  a  unit  dipole  when  on  the 

boundary  between  air  (Region  0,  z  ±  0)  and  a  dielectric  or  conducting 

half-space  (Region  1,  z  2  0)  is  well  known.  The  radial  component  of 

the  electric  field  (f ,  <$  ,  0)  includes  a  wave  of  the  type  known  as 

a  lateral  wave  which  travels  in  the  air  along  the  boundary  and  a 

direct  wave  that  travels  in  the  adjacent  half-space.  The  lateral  wave 
-  ifto?  . 

in  Region  0  has  the  coefficient  e  with  associated  amplitude  factors, 

lkl? 

the  direct  wave  in  Region  1  has  the  coefficient  e  with  associated 
amplitude  factors.  Because  the  wavenumber  k  differs,  these  two  waves 
propagate  at  different  velocities.  Unfortunately,  these  known  formulas 
for  the  air-dielectric  interface  are  not  valid  for  a  three-layered 
region,  like  an  integrated  circuit. 

Corresponding  new  formulas  were  derived  for  the  horizontal  electric 
dipole  on  the  boundary  of  a  three-layered  region  shown  in  Fig.  8.1. 

The  mathematical  details  are  too  lengthy  to  be  presented  in  this  report 
and  will  be  published  elsewhere.  Only  the  major  research  results  and 
their  relevance  to  the  coupling/interference  problem  in  very-high-speed 
integrated  circuits  are  discussed  here. 

The  following  practical  approximations  were  used  to  derive  the 
new  formulas  for  the  fields: 

|ko2|<<  |kx2|<<  |k22J,  jkjj?!  c  1,  (8.1) 

where  the  wavenumber, 

1/2 

kj  =  +  ie<j  =cdC +  *  (8.2) 

and  2nf,  the  radian  frequency.  In  (8.2)  the  quantities €, 
and  a-  have  their  usual  meanings.  In  (8.1),  |k^$  |  is  interpreted  as 
the  electrical  thickness  of  the  dielectric  layer  whose  physical 
thickness  is  &  . 
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The  general  form  for  the  radial  electric  field  at  the  air-dielectric 
surface  due  to  a  horizontal  electric  dipole  (see  Fig.  8.1)  is. 
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This  includes  contributions  from  a  surface  wave  that  travels  in  the 

ik-.  P 

dielectric  region  1  with  the  propagation  factor  e  v  and  a  lateral  wave 

that  travels  along  the  surface  in  the  air  with  the  propagation  factor 
iko0 

e  '.  These  two  waves  have  different  velocities  of  propagation. 

The  quasi-static  or  near  field  part  of  (8.3)  is. 
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(8.4) 


At  large  distances,  the  lateral  wave  dominates  and. 
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The  general  form  for  the  phi  component  of  the  electric  field  at 
the  air-dielectric  interface  due  to  a  horizontal  electric  dipole 
(see  Fig.  8.1)  is, 
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The  quasi-static  or  near-field  part  of  (8.6)  is, 

^o  sin(p  _ 

El^(^,cp,  0)~ - [1  +  2ik^],  |k^p  <-<  1. 

2?rk12  p3 

At  large  distances  the  far  field  has  the  form. 
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,  0)  oj  -  sin<|)  - - [1  -  kJ],  |k1p  I  <<  1. 
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Although  the  general  formulas  for  the  components  of  the  electric 
field,  (8.3)  and  (8.6),  are  complicated,  some  important  observations 
relevant  to  coupling  and  interference  in  integrated  circuits  can  be  made. 
First,  reducing  electric  field  strength  in  the  region  near  the  current 
element  is  desireable  because  coupling  and  interference  are  proportionally 
reduced.  The  near  field  terms  in  (8.4)  and  (8.7)  each  have  the 
coefficient, 


'l  [///^  +  iO^/w)]172 


(8.9) 


Therefore,  increasing  the  permittivity  e  of  the  dielectric  layer 
reduces  field  strength  (and,  consequently,  coupling  and  interference 
with  nearby  conductors  and  circuits).  At  high  frequencies  (J^/co  becomes 
smaller  tending  to  increase  field  strength.  Increasing  the  conductivity 
of  the  dielectric  layer  to  compensate  this  may  be  a  poor  choice  because 
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signal  losses  in  the  dielectric  will  increase  and  heating  could  be 
problematic. 

At  larger  distances  away  from  the  current  element  the  lateral  wave 
has  the  coefficient. 
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(8.10) 


Again,  field  strength  is  reduced  by  increasing  the  dielectric's 

permittivity  G Coupling  and  interference  also  can  be  reduced  at 

larger  distances  from  the  current  element  by  reducing  the  physical 

thickness  X  of  the  dielectric  layer  as  can  be  seen  from  (8.10).  Because 
2 

kQ  appears  in  the  numerator  of  (8.10),  increasing  kQ  would  proportionally 

increase  the  strength  of  the  lateral  wave  and  coupling  or  interference 

due  to  it.  This  means  that  using  a  material  with  G  >£  in  place  of 

o 

air  in  Region  1  is  not  good  from  the  standpoint  of  interference/coupling 
due  to  the  lateral  wave. 

The  general  expression  for  the  radial  electric  field  contains 
several  l/k2  terms  which  can  be  made  smaller  by  increasing  the  con¬ 
ductivity  0~2  of  the  ground  plane  as  can  be  inferred  from  the  relation, 

k2  =W[A(€2  +  ia’2/°3^l/2’  (8.11) 

In  summary,  from  the  field-theoretic  path  coupling  model  based  on 
a  study  of  fields  generated  by  a  horizontal  current  moment  parallel  to 
the  air-dielectric  surface  of  the  three-layered  problem,  coupling  and 
interference  generally  can  be  minimized  by:  (1)  increasing  the 
permittivity  €1  £f  the  dielectric  layer;  (2)  decreasing  the  thickness  £ 
of  the  dielectric  layer;  (3)  increasing  the  conductivity  of  the  ground 
plane  below  the  dielectric  layer.  These  results  are  consistent  with 
that  from  the  previous  section,  namely,  increase  path  capacitance  to 
minimize  coupling  and  interference. 

The  interesting  new  research  result  in  this  section  is  the 
existence  of  two  mechanisms  of  coupling  from  a  current  moment  on  the 
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dielectric  surface:  one  due  to  a  lateral  wave  and  the  other  due  to 
a  surface  wave.  Because  the  velocity  of  propagation  of  one  wave  is 
different  from  the  other,  these  interfering  signals  arrive  at  some 
point  on  the  surface  of  the  dielectric  at  different  times.  This  may 
be  critical  at  very  high  clocking  rates. 
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